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SUMMARY 

A previous analysis of turbulent heat transfer and flow with vari- 
able fluid properties in smooth passages is extended to flow over a flat 
plate at high Mach numbers . Velocity and temperature distributions are 
calculated for a boundary layer in which the effects of both frictional 
heating and external heat transfer are appreciable. The viscosity and 
thermal conductivity are assumed to vary as a power of the temperature, 
while the Prandtl number and specific heat are taken as constant. Skin- 
friction and heat-transfer coefficients are calculated and compared with 
the incompressible values. The relation between boundary- layer thickness 
and distance along the plate is obtained for various Mach numbers. The 
analytical results are compared with representative experimental data. 


INTRODUCTION 

Okie current emphasis on high-speed flight has caused much interest 
in research on compressible boundary layers . The skin friction in high 
Mach number flight constitutes a large part of the total drag. There- 
fore, the accurate prediction of skin friction is desirable for the de- 
sign of high-speed aircraft. Prediction of heat-transfer coefficients 
in high Mach n umb er flow is also important, because frictional heating 
of the surface necessitates cooling to prevent structural failures. 

The prediction of laminar boundary layers from the basic equations 
of momentum, energy, and continuity has reached a high state of develop- 
ment. A considerable amount of analytical work an turbulent boundary 
layers has al so been carried out. In the turbulent case, however, the 
results of the various analyses disagree markedly because of the differ- 
ent assumptions made by the various authors. These a n a l yses are reviewed 
in references 1 to 3. The introduction of assumptions into the treatment 
of turbulent bo undar y layers is at present unavoidable, since solving the 
problem from the instantaneous equations of momentum, energy, and conti- 
nuity alone is not yet possible. In seme respects, however, the model 
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used for solving the problem might be improved. In nearly all the anal- 
yses , the flow is divided into a laminar region, "where turbulence is sup- 
posed to be absent, and a fully turbulent region. The effect of varia- 
tion of fluid properties on the laminar region is generally neglected. 
Measurements of turbulent velocity profiles indicate that considerable 
turbulent shear exists within the so-called laminar layer (ref. 4), so 
that a more realistic model for the region close to the wall than that 
used in previous analyses is desirable. 

A somewhat improved treatment of the region close to the wall is 
given in references 4 to 6, where the effects of turbulence and of vari- 
able fluid properties in this region are considered. In the region away 
from the wall the von Karmen similarity expression has been considered 
the most reasonable expression available (ref. 7). In reference 8, fully 
developed turbulent flow and heat transfer in smooth passages for air 
with variable properties are analyzed, and the results agree well with 
experimental data. The analysis is extended to the entrance regions of 
passages and to high Prandtl numbers In references 9 and 10, where good 
agreement with experiment is again obtained. Since the analyses apply 
well to entrance regions, the assumptions made in the analyses should 
apply also to a compressible boundary layer. ""The analysis is extended 
to flow and heat transfer in a boundary layer at high Mach numbers in 
this paper. (Seme preliminary results were presented in ref. 11.) The 
variation of properties due both to frictional heating and to external 
heat transfer is considered. The viscosity and thermal conductivity are 
assumed to vary as a power of the temperature, while the Prandtl number 
and specific heat are taken as constant . 
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constant 

ratio of diffusivities, e - n / e 

constant 

constant 

friction coefficient, 2 t v / pgUg 

specific heat of fluid at constant pressure 

constant 

exponent for viscosity variation with temperature, 
for air 
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taken as 0.68 
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constant 

enthalpy 

heat-transfer coefficient, q w /(t w - t ) 
thermal conductivity 

Mach number based on free -stream' properties and velocity. 


WrRtJ 
constant, 0.109 
Prandtl number, c p^/k 

heat transfer In y-direction per unit time per unit area 

perfect gas constant 

Reynolds number based on x, xugpg/(ig 

Reynolds number based on 9 , SUgpg/dg 

Stanton number, h/c^UgPg 

total temperature, t + (u^/2c } , deg abs 


total-temperature parameter. 


K - Vv 1 - wy) 


static temperature, deg abs 

(t w - t)CpT w 

temperature parameter, ■ — = 

^wV T w/ p v 

2 (V " t ) c p p v 

temperature parameter, — 

v 

velocity in x- direct ion 


1 - (t/fcg) 

P 

1 - (Vjy) 

a 


velocity parameter, u/ ->J t^/ p v 
velocity in y-direction 
longitudinal distance along plate 



4 


NACA M 4262 


distance perpendicular frcm plate 


wall distance parameter. 




lowest value of y + for which equation for region away frcm wall 
applies 

frictional -heating parameter, T w /2c^t w p w . 

^wV^/pw 

heat-transfer parameter, ~ 

c p t w T w 

ratio of specific heats, taken as 1.400 for air 


flow boundary-layer thickness 


flow boundary-layer-thickness parameter. 




thermal boundary-layer thickness 


thermal boundary-layer-thickness parameter. 


eddy diffusivity of momentum 




^Sf/Pw 


eddy diffusivity of heat 


temperature-recovery factor, — g 


"^av " *8 




momentum thickness, / -2 — — ( 1 — — ) dy 

4 ) p 5 U S \ 


momentum- thickness parameter, 


e V T w/Pw 


constant, 0.36 

viscosity 

density 

shear stress, force per unit area 
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Subscripts: 

aw pertaining to adiabatic wall conditions 

i incompressible; constant fluid properties 

w pertaining to wall 

6 pertaining to edge of boundary layer or free stream 

1 pertaining to edge of wall layer 

Superscripts : 

* reference 

1 pertaining to fluctuations from. time average except in t + ' 

time- averaged value 


ANALYSIS AND DISCUSSION 


Basic Equations 

The instantaneous velocities, temperature, and fluid properties in 
the equations of momentum, energy, and continuity can be divided into 
mean and fluctuating components . If time averages are taken, the follow 
ing equations for shear stress and heat transfer, applicable to flow in 
a boundary layer, are obtained (appendix A): 


r = H ^ - P u'v' (l) 

q = “ k ly + pC p t,v ’ " ^ ly + ^ u ' v ’ ( 2 ) 

where constant specific heat is assumed. The bars denote time averages, 
and the primes indicate fluctuating components. Equations (l) and (2) 
are the same as equations (A9) and (A14) in appendix A if the bars over 
the time-averaged velocities, temperatures, and properties are dropped. 
The various terms in equations (l) and (2) may be interpreted as follows 



molecular shear stress 


-p u'v 1 turbulent shear stress 



molecular heat transfer 
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pCp t'v* turbulent heat transfer 
du 

-up. molecular dissipation 

up u'v' turbulent dissipation 

Equations (l) and (2) suggest the form of the turbulent transfer 
equations but contain the unknown quantities u'v’ and t'v', so that 
assumptions must be made before solutions can be obtained. For making 
these assumptions it is convenient to introduce the relations 


u'v' = -€ -=~ and t 'v 1 = -e*. ~ 

<3y h ay 

where e and are the eddy diffusivities for momentum and heat trans- 

fer, the values of which depend upon the amount and kind of turbul e nt 
mixing at a point. When these relations are introduced, equations (l) 
and (2) become 

T “ (t* + P £ ) ^ (3) 


q = - (k + pc p e h ) ^ 


The physical significance of e and lies in the fact that e/(ii/p) 

is the ratio of turbulent to molecular shear stress (ref. 12), and 
e-h/(k/pc ) is the ratio of turbulent to molecular heat transfer. Equa- 
tions (3) and (4) can be written In dimensionless form as 



and 

= 1 + _P. a _4_\dt^ _ 2 _; u +fjL + £4_W 

p w v^jdy+ P \V^ Pw p w/ p w ; dy+ 


( 6 ) 


The subscripts w refer to values at y = 0$ that is, at the wall. The 
quantity a is a frictional-heating parameter that is an indication of 
the variation of properties due to frictional heating, and j3 is a heat- 
flux parameter that is an indication of the variation of properties due 
to heat transfer. The parameter a is always positive or zero, a value 
of zero characterizing low-speed flow (i.e., Mg = 0) . A zero value of 
6 refers to a vanishingly small heat transfer or an insulated plate. 
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if 


A positive value of p indicates heating of the fluid, while. negative 
p means that the fluid is being cooled. It is sometimes convenient to 
write equation (6) in the following alternative dimensionless foim: 


— — = 1 P a e \ dt + _ + fji_ _p_ e \ du + 

a q w \V Pr w Pw * <3y + Pw tV/ p wJ dy+ 


(7) 


This equation is particularly convenient when p = 0, for which case 
equation (6) becomes indeterminate . 


Expressions for Eddy Diffusivity 

In order to make practical use of equations (5) to (7), the eddy 
diffusivity e must be evaluated for each portion of the flow. For this 
purpose the boundary layer is divided into two portions termed the "re- 
gion away from the wall" and the "region close to the wall." 

Region away from wall . - In the region away from, the wall, it is 
assumed that the turbulence at a point Is a function mainly of local con- 
ditions - that is, of the relative velocities in the vicinity of the 
point (ref. 13). This is probably not a good assumption in the region 
near the edge of the boundary layer, where considerable diffusion of the 
turbulence occurs (ref. 14) and, in addition, the turbulence is inter- 
mittent. However, in that outer region the velocity or temperature 
gradients are so small with respect to these gradients nearer the wall 
that the error in calculated velocities or temperatures should not be 
large. A Taylor series expans ion for u as a function of transverse 
distance, then, indicates that e is a function of du/d y, d 2 u/dy 2 , 
d^u/d y^, and so forth. If, as a first approximation, e is considered 
as a function only of the first and second derivatives, and dimensional 
analysis is applied. 



This expression was obtained by von K&rman and is generally known as the 
K&rm£n similarity hypothesis (ref. 7). The constant x is to be deter- 
mined experimentally. 

Region close to wad 1 . - In the region close to the wall it is as- 
sumed that e is a function only of quantities measured relative to the 


* 
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wall - that is, of u and y. This assumption includes, to a first 
approximation, an effect of the derivative du/dy. Since the flow be- 
comes very nearly laminar as the wall is approached, the first deriva- 
tive approaches the value u/y and hence may be omitted, since u and 
y already appear in the functional relation. By using d ime nsional 
ana-lysis. 


e = e(u,y) = n 2 uy (9) 

where n is an experimental constant. 

Equations (8) and (9) can be considered as reasonable first approx- 
imations for e. Whether these approximations are adequate or not can 
at present be determined only by experiment. 

Determination of experimental constants . - The constants n and x 
were determined from pipe data in which the properties were essentially 
constant. Equation (5), with equation (8) or O), was integrated (con- 
stant properties and t) for the regions close to and away from, the wall 
in reference 4. The molecular shear stress was neglected in the region 
away from the wall, and the well-known Karman-Prandtl logarithmic equa- 
tion was obtained in that region. In matching the two solutions it was 
assumed that the velocity is continuous at the junction of the two 
regions . 

The integrated equations (ref. 4) for the regions close to and away 
from the wall are platted in figure 1 with the constants n = 0.109 and 
x = 0.36 determined from pipe data (refs. 4 and 14). The data indicate 
that the equation for the region close to the wall applies for y+ < 26, 
and the equation for the region away from the wall applies for y + >26. 
Included in the plat are data for a low- speed boundary layer with zero 
pressure gradient from reference 15. The agreement with the curve is 
satisfactory . 

The values for the constants n = 0.109 and x = 0.36 should apply 
to flow with variable as well as constant properties if the basic assump- 
tions made for e in the preceding sections apply to variable properties; 
that is, if e = e(u,y) close to the wall and e = e(du/dy, d^u/dy^) 
away from the wall. The constant y+, however, requires further consid- 
eration and is discussed in the next section. 


#• 


£ 

a 


as 


^Reference 10 shows that the kinematic viscosity has an effect on 
e in the region very close to the wall. However, that effect becomes 
important only for heat or mass transfer at Prandtl or Schmidt numbers 
appreciably greater t han 1. 
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Additional Assumptions 

In addition to the assumptions for eddy diffusivity discussed in 
the preceding section, several additional assumptions must he made for 
solving equations (5) to (7). 

Variation of properties -with temperature . - For gases, the viscosity 

varies approximately as t^, where d has an average value of 0.68 for 
temperatures between 0° and 2000° F. The Prandtl number (Pr = 0.73) and 
specific heat Cp are assumed constant, because their variations with 

temperature are of a lower order of magnitude than the variations of the 
other properties. If and Prandtl number are considered constant, 

the thermal conductivity k will vary with temperature in the same way 
as the viscosity, or as t* 1 . For constant pressure across the boundary 
layer, the density p is inversely proportional to t. 

With the preceding assumptions, the property ratios in equations 
(5) to (7) can be written as 

JL 

^w 


From the definitions of |3 


or, if equation (7) rather 


The property ratios in equations (5) to (7) can therefore be written in 
terms of p and t + or a and t + . 

Variations of t and q across boundary layer . - The momentum 
equation (A7) indicates that, for a flat plate (zero pressure gradient), 
dx/dy = d(pdu/dy) /dy = 0 at the wall. Since r is zero at the 
edge of the bo undar y layer, the actual variation of t across the 
boundary layer might be expected, in general, to lie between a linear 
variation = 1 - (y/8) ) and t/t^ = 1. Data on low-speed isothermal 



and t + , 

5T- = 1 - |3t + (12) 

z v 

than equation (6) is used, 

— = 1 - at + ' (13) 
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flow over a flat plate (ref. 15) show that this type of variation does 
exist, except in a narrow region near the edge of the boundary layer. 

For determining the sensitivity of the velocity or temperature profile 
to shear-stress variation, it should therefore he sufficient to compare 
the profiles for a constant and for a linearly varying shear stress. 
Appendix B shows that x/ t w = q/q w for a flat plate if the Prandtl num- 
ber is 1. 

Figure 2 shows u + or l 4- plotted against y + for a Prandtl num- 
ber of 1 for both a constant and a linearly varying shear stress and 
heat transfer, where T 4 is the total-temperature parameter. Curves 
are shown for 6 = 0 and a - 0, 0.003, and 0.008, which cover much of 
the range of Mach number and Reynolds number of interest. The equations 
for calculating the curves are given in appendix D. The equation for 
the region away from the wall was taken to apply for y* > 30 rather 
than > 26 when the shear stress was variable, in order to give better 
agreement with the data for constant properties. The curve for a = 0.008 
is cut off at the point shown because the Mach number becomes infinite, 
as can be seen from equation (D3) (for au +2 = l) . The curves indicate 
that variable shear stress and heat transfer have but a slight effect on 
the velocity and temperature profiles. Similar curves were obtained in 
figure 11 of reference 6 for p ^ 0 and a = 0. The same conclusions 
should apply to Prandtl numbers differing slightly from 1, so that the 
effects of the variations of x and q across the boundary layer are 
neglected for solving equations (5) to (7). 

Ranges of applicability of equations for flow cl ose to and away 
from wall . - It was determined from, the data for constant properties 
that the lowest value of y + for which the equation for the region 
away from the wall applies is y * = 26 when the variation of shear stress 

with y is neglected and the molecular shear stress is neglected in the 
region away from the wall. The question arises as to how y^ varies 

when the properties are variable. The simplest assumption is that yj 

is constant and equal to 26. This assumption, which implies that the 
wall properties govern the thickness of the wall layer 

(y^ 2 yi Vv^/ (ll^/p ^) ) > is similar to vcn Karmen's assumption (ref. 

16). Figure 12 of reference 6 shows that essentially the same curves 
are obtained when the molecular shear stress is neglected in the region 
away from the wall as when it is considered, the difference being that, 
when the molecular shear stress is included, y* has the constant value 
of 16 rather than 26. 

Another assumption, which might be somewhat more reasonable tban 
assuming y* constant, is that y£ occurs at a given constant ratio 

of turbulent to molecular shear stress e/(p/p). That is, the turbu- 
lence chan ges from that described by equation (9) to that described by 
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equation (8) when the ratio of turbulent to molecular shear stress 
reaches a certain value. In this case the more complete equations are 
used for the region away from the wall, in which the molecular shear- 
stress and heat-transfer terms are retained and the slopes of the equa- 
tions for flow close to and away from the wall are matched at y£ 

(y* = 16 for p = a, = 0) . 

In figure 3, u + or T 4 " is plotted against y + for a Prandtl num- 
ber of 1 using the two assumptions for y£ discussed in the preceding 
paragraphs. Curves are shown for p = 0 and a, = 0, 0.003, and 0.008. 

The equations for calculating the curves are given in appendix C. The 
curves indicate that the velocity and temperature profiles are apparently 
insensitive to the assumption used for y^. Similar results were ob- 
tained in figure 13 of reference 6 for P ^ 0 and a = 0. The simpler 
procedure of neglecting the molecular shear stress and heat transfer in 
the region away from the wall and assuming y* = constant => 26 is there- 
fore adopted in the following calculations. 

Ratio of eddy diffusivities for heat and momentum, transfer . - In 
most analyses the ratio of eddy diffusivities a that occurs in equa- 
tions (6) and (7) is set equal to 1; that assumption has given heat- 
transfer coefficients in good agreement with experiment (ref. 8). It is 
* of interest that Prandtl 's mixing-length theory, which assumes that a 

turbulent particle moves a given distance and then suddenly mixes with 
the fluid and transfers its heat and momentum, gives a value of a = 1. 
Although the actual turbulence mechanism may be more complicated than 
indicated by that theory, it does indicate that a value of a on the 
order of 1 is not unreasonable . 

In the present analysis the assumption of a = 1 is retained, but 
in some cases the calculations are also carried out for a = 1.07 in 
order to determ in e the effect of varying a. A ratio of diffusivities 
of 1.07 was obtained from seme preliminary experiments on recovery fac- 
tors for f ully developed flow in a tube. 


Velocity and Temperature Distributions in Boundary Layers 


For obta inin g velocity and temperature distributions close to the 
wall, equations (9) to (13) are substituted into equations (5) to (7). 
Equations (5) and (6) become, in integral form, with t/t w = q/q w = 1, 



3y^ 

(1 - pt + ) d + - — - n 2 u + y + 

1 - pt + 


( 14 ) 
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Equations (14) and (l5) can be solved simultaneously by iteration^ that 
is, assumed relations between u + and y + and t + and y+ are sub- 
stituted into the right sides of the equations, and new values of u + 
and t + are calculated by numerical integration. These new values are 
then substituted into the right sides of the equations and the process 
is repeated until the values of u + and t + do not change appreciably. 
Equations (14) and (15) give the relations among u + , t + , and y + for 
various values of a and 6 for flow close to the wall (y + < 26). For 
|3 = 0 and a =4 0, t + becomes infinite, so that equations (5) and (7) 
must be used. These equations, with equation (13), become 



Equations (16) and (17) are solved similarly to equations (14) and (15). 


In the region away from the wall, the molecular shear stress and 
molecular heat transfer are neglected. Dividing equation (6) by equa- 
tion (5) gives, with t/t^. = q/q w = 1, 

1 + 2 - u T = a — - (18) 

P du + 

Integrating equation (18) from y* to y + gives 





a 

pa 




( 19 ) 
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> i . 


From equations (ll), (l2), and (l9). 


1 - stt + ^ + 2 „+2 . . 2 u +2 

J- a a 


( 20 ) 


a a 


Substitution of equations (8) and (20) into (5) and one integration give, 
for the region away from the wall, 


v <iu 

Ka T-+ = e 
dy^ 


sin 


a/ a 


■1 2ou + +p H 

-J p 2 +4a( a- apt^+gu^+au^ 2 )J 


By letting 


x _i 

z = — : sin 


a/ a 


+ 

2au+p 


+ 4a(a - apt^ + pu£ + <xuj^) 

and integrating equation (2l), 


+ ^b S ( -JoJe. , -Ja/a . i/a/a \ 

f = 77 JL \ z Z J 


+ K 


ax" 


( 21 ) 


( 22 ) 


(23) 


The constant K is evaluated in the usual way by letting 
du + /dy + = °° at y + = 0 in equation (2l) (ref. 7) and substituting (2l) 
into (22) and (23) at y + = 0.^ By using this procedure, K = 0. To 
determine K^, set u + = u^ when y + = y£. Then, 



where z is given by equation (22) and z^_ is the value of z at 
y£ = 26. Equations (22) and (24) give the relation between u + and y + 


p 

This assumption can be avoided by including the molecular shear 
stress and heat transfer in the region away from, the wall and evaluating 
K by assuming a continuous velocity derivative at y* (fig. 12, ref. 6). 
This assumption gives essentially the same results as that made in the 
text . 
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for various values of a and (3. The quantity t + can then he calcu- 
lated from equation (19) . 


For 6 = 0 and a ^ 0, t + becomes infinite, and t + ' = (p/a)t + 
must he used. Equation (19) becomes, in terms of t + , 


t + ’ = <’ - i < 
1 aa l 


u l 2 B 4 - u +2 

-±- + u + + 
a aa a 


(25) 

Equations (22) and (24) apply to the case for 6 = 0 and a ^ 0 if at* 
is substituted for pt£ in equation (22). 

For a = 0, equation (23) hecomes indeterminate, and equation (2l) 
for zero frictional heating from reference 6 can he used. 


Typical velocity and temperature distributions for various values 
of the frictional-heating parameter a and of the heat-flux parameter 
p are presented in figures 4 to 6. Positive values of p correspond 
to heat addition to the air; negative values, to heat extraction. The 
curves of u + against y + (fig. 4) indicate considerable flattening of 
the velocity profile as either a or p increases positively. This is 
caused by the decreasing temperatures in the outer regions of the bound- 
ary layer compared with the wall temperature when either the Mach number 
is high (high a) or the heat transfer from the surface to the air is 
high. Thus, the density is higher in the outer regions of the boundary 
layer, with consequent flattening of the profile (eq. (5)). Negative 
values of p produce the opposite effect. For certain combinations of 
a and 6 (with p negative), the effect of 6 on the curves should 
tend to cancel the effect of a, and the resulting profile should not 
differ greatly from the a = p = 0 curve. The curve for a = 0.002, 

P = -0.05 in figure 4 is close to the curve for a = p = 0. Included 
in figure 4 for comparison are experimental data from reference 17 for 
an a of 0.00176, p = 0, and a corresponding Mach number Mq of 2.82. 

The data are in reasonable agreement with the predicted profiles. 

J_ j. 

In figure 5 T is plotted against y^ for various values of a 
and p . The total- temperature parameter T* is plotted rather than t + , 
because the trends are seme what more consistent, although some crossing 
over of the curves occurs even with T + . For calculation purposes, a 
better representation can be obtained by plotting T 4 " against u + . The 
quantity T* is related to t + by the relation 


T 1 " = t + - # u + 2 


(26) 


QTP^ 
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Skin-Friction Coefficients 
The skin-friction coefficient is defined as 


2t. 


= 


w 


Ps u i 


(27) 


where the subscript 8 refers to values outside the houndary layer. 
Equation (27) becomes, in dimensionless form.. 


p w 2 


2(1-0*$') 2(1 - pt$) 




+2 

*8 


(28) 


For comparison with experimental data, it is convenient to introduce the 
momentum, thickness. 


which in dimensionless form is 


e = 

/ n 

— f 1 - 




V. 

L p s 

“s \ 

"8 J 



form 

is 





, + '< 
atg , 

. r 

1 

u + 

A 

_ u+ 

X 

1 - at + 


K 

+ 

U 5. 


(29) 


dy* 


(30) 


Then the Reynolds number based on the momentum thickness and free-stream 
properties is 


„ _ eU 8 P 8 n + + ^w P 8 

Re Q = = Q u~ 

0 p S ^ p 8 p w 

where the property ratios are obtained frcm equations (10) to (13), 
Mach number for a perfect gas is 


% “ 


*8 



If values are given to a, f3, and 


2a 

(r - i;(t s /t w ) 


(31) 


The 


(32) 


8 + , where 


8 + is the value of 


y + at the edge of the flow boundary layer, then values of ug, tg, and 
so forth can be read frcm curves similar to those in figures 4 to 6. 

Values of C f , Re 0 , and Mg can then be calculated frcm. equations (28), 

(3l), and (32). This procedure assumes that the thermal and flow boundary 
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layers are of equal thickness . From the calculations in a later section, 
where relations between boundary-layer thickness and distance along the 
plate are calculated, it can be shown that this is a good assumption for 
gases when the thermal, and flow boundary layers begin at the same point. 
For the case of Pr = a = 1 the assumption holds exactly, as can be seen 
by substituting u + = T* into equations (42) and (43), which are then 
identical . 

Predicted skin-friction coefficients are plotted against Reg in 
figure 7 for various values of Mach number for an insulated plate 
(0=0). These curves are for e^/e = a = 1. The effect on the curves 

of changing a to 1.07 was negligible. The values of Ch, decrease 
considerably as Mach number increases. Included in the plot are experi- 
mental data of a number of investigators for Mach numbers up to 4.93. 

In general, the data are in good agreement with the predicted curves . 

The ratio of the friction coefficient to the inc expressible coeffi- 
cient is plotted against Mach number for various values of Reg for 

0 = 0 in figure 8. The values of Vf/Vf } ± decrease with Re 0 , but at 

a decreasing rate. For comparison purposes the analytical curve of 
Cf/Cf against Mg for a value of Re 0 of 6000 is plotted in figure 
9 together with data taken near this value of Re . 

C7 

If heat transfer occurs between the plate and the stream, it iB 
convenient to specify the ratio of the actual wall temperature to the 
adiabatic wall temperature for a given Mach number and Reg. For an in- 
sulated plate the adiabatic wall temperature may be written as 

2 

*«• " *8 + 1 I" < 33a > 

where r\ is the temperature -recovery factor, the calculation of which 
is discussed in the next section. Equation (33a) can be written in 
dimensionless form as 



1 - ptj + 


(33b) 


Figure 10 is similar to figure 7, except that the plate is now 
cooled (t w /t aw = 0.5). The trends are similar to those of figure 7, but 

all the curves are displaced upward. This increase in friction coeffi- 
cient was also obtained for flow in a tube with cooling (ref. 5). Also 
included in this figure are wind-tunnel data (ref. 18) obtained at high 
Mach number using nitrogen as the working fluid. The agreement with 
theory appears to be within experimental error. 


enro^ 
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Figure 11 is similar to figure 8 but is for a value of 0.5 for 
t v /t aw . It is of interest to note from both figures 8 and 11 that the 

percentage effect of varying Reynolds number is much greater for the 
higher Mach numbers . Thus,, figure 8 indicates that for a Mach number of 
20 the value of C^/Cf for Re^ of 10® is less than half that for 

Re 0 of 10®. 


Reference Temperatures and Extension of Results to 
Greater Cooling Rates 

The defining equation for reference temperature is 

t* = tg + C(t w - tg) + D(t av - t w ) (34a) 


where C and D are constants to be evaluated from theoretical or ex- 
perimental results. Dividing equation (34a) by tg and assuming that 
the recovery factor is constant at 0.88 and that x — 1-40 result in a 
more useful form of the equation: 

t * 

= (1 - C) + 

t 8 

According to the concept of reference temperature, the relation be- 
tween incompressible friction factor and Reynolds number should hold for 
variable-property flow if the properties are evaluated at the reference 
temperature . The results of this analysis could not be represented ac- 
curately in such a manner. It was necessary to write the incompressible 
relation in the form 

C - E - 
*>*■ Reg-0744 


(C - D) + D 


aw 


1 + 0.176 


«§) 


where E is a constant. It should be noted that this is not the true 
incompressible relation and 1 b used only for reference temperature pur- 
poses . If the properties are evaluated at the reference temperature, 
the result is 


C 


f 


E 

R e 0 • 0744 


'^_*\-0.875 

ft) 
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Dividing this equation by the previous one gives 


C 


f 




(34b) 


For Reg of 10^, evaluating the constants C and D in the refer- 
ence temperature equation from the results in figures 8 and 11 and equa- 
tion (34b) gives C = 0.56 and D = 0.184. The reference temperature 
can then be written 


= 0.44 + 


^0.376 — - + 0.184^1 + 0.176 m|) 


(34c) 


The results of the use of equations (34b) and (34c) are shown as dashed 
lines in figures 8 and 11. Thus, by use of equations (34b) and (34c) it 
should be possible to extend the results of this analysis to values of 
t M / t aw other than 1.0 and 0.5 if the value of Reg is near 10^. 


An estimate for lower Reynolds numbers may be obtained by first 
using the preceding method to find the friction factor at Re 0 = 10^ 
and then finding the ratios of friction factor at the desired Reynolds 
number to that at Reg = 10 5 from figures 8 and 11. These ratios can 


then be interpolated or extrapolated to the desired value of t w /t aw . 
This procedure can be justified since the ratio does not vary greatly 
with t w /t aw . 


Stanton Numbers and Recovery Factors 

The Stanton number based on the difference between the wall and the 
adiabatic wall temperature, with properties evaluated at the free-stream 
temperature, is given by 


St = 


h 

c p u SPS 


1 


. + + PS 




- r)ccug 


(35) 


where equation (33a) is used. The temperature and density ratios are 
determined from equations (ll) and (12) or (13) . For 6 = 0 and a, / 0, 
equation (35) becomes indeterminate. For that case set 

1 - (tg/t w ) = atg' and t* = (<x/6)t* . With these substitutions equa- 
tion (35) becomes 
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A 


St = 


Fran, equations (7) and (5) (q/q. w = t/t v - l). 


(36) 


t + ’ - P 
*8 “a 


<3y 


2u + dy + 


■4 \ Pr w P w P'w/p-w ^0 \ Pr v P W P\/ P w 


(37) 


The second integral in this equation can be replaced by (tg*) But 

equation (33a) can be written in dimensionless form as 


^ 2 - (t sVo w 

Substituting equations (37) and ^38) into (36), with the second integral 
in equation (37) replaced by (t* )p = o> gives 


St 


1 



dy + 

Jl _i_ + _e_ 


a 


Pw/ p w 


(39) 


For evaluating equation (39) in the region close to the wall, 
e/([i w /p v ) = n 2 u + y + . For the region away from, the wall, e could be ob- 
tained from, equation (8) . However, it is more convenient to obtain e 
from equation (5), which for the region away from the wall becomes 

p e 1 

p w Py;-/ p w du + /dy + 

Equation (39) can be used for 6 = 0 or p / 0. For given values of 
8 + , a, and 0, values of Stanton number, Mach number, t w /t aw , and Reg 

can be calculated from equations (39), (32), (33b), (30), (3l), and (10) 
to (13). 

Predicted Stanton numbers are plotted against Re^ for various 
Mach numbers for t w /t aw = 1 in figure 12. The case of t w /t aw = 1 is 

a limiting case that can be approached as closely as desired by making 
the heat flux small. When t w /t av = 1, there is no effect of variable 
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properties due to beat flux. The Stanton numbers in figure 12 show- 
trends similar to those of the friction coefficients in figure 7. In- 
cluded in the figure are experimental data for low heat flux obtained by 
a number of investigators. In general, the data are in good agreement 
with the predicted curves. The curves in figure 12 are for e^/e = a *> 1. 

Similar curves for an a of 1.07 were 3 to 5 percent higher for a Mach 
number of 0, but the difference decreased at higher Mach numbers . The 
curves for a = 1 are in slightly better agreement with the data than 
those for a = 1.07. 

The ratio of Stanton number to the incompressible Stanton number is 
plotted against Mach number for various values of Re 0 for t^/t^ = 1 

in figure 13. These curves are very nearly the same as those for 
C^/Cf i in figure 8. 

Figure 14 is similar to figure 12, except that t /t =0.5. As 

was the case for the friction coefficients in figure 10, the Stanton 
numbers increase as t w /t aw decreases . The corresponding plot of 

St/St^ against Mach number for various values of Reg and t w /t = 0.5 
is shown in figure 15. * 

Temperature-recovery factors, as calculated from equation (38), are 
shown in figure 16 for Mach numbers from 0 to 8. Curves are shown for “ 

a = 1 and a = 1.07. The curves for a = 1.07 are in somewhat better 
agreement with most of the experimental data than those for a = 1. 

This does not mean that an a of 1.07 should be used for calculating 
heat transfer or Stanton numbers. According to Reichardt’s hypothesis, 
the value of a should be close to 1 at the w all and increase with dis- 
tance from the wall (ref. 12). The temperature profiles for calculating 
heat-transfer coefficients are very steep near the wall, so that the im- 
portant charges with distance take place near the wall where a is close 
to 1. In the case of recovery factors, however, the plate is insulated, 
so that the temperature gradient is zero at the wall. The gradients near 
the wall will therefore be smaller than in the case of heat transfer, and 
imp ortant changes of temperature with distance might occur in regions 
away from the wall where a is somewhat greater than 1. 

Figures 17(a) and (b) show the curves of the Reynolds analogy factor 
2St/C f against Mach number for various Re 0 , for t^/t^ equal to 1.0 

and 0.5, respectively. If Reynolds’ analogy held strictly (Pr = a = 1), 
the Stanton number would be equal to cue-half the friction factor and 
2St/c^ would be unity. Figure 17 shows a variation of the Reynolds 
analogy factor over the range 1.065 to 1.280. In general, 2St/C|. in- 
creases with Increasing Mach number and with decreasing Reynolds number 
and increases slightly with decreasing t w /t aw at the higher Mach 
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numbers . These results are in approximate agreement with those of 
Rubesin (ref. 19), who estimated that 2St/C.p would he in the range 
1.18 to 1.21 at least up to Mg of 5. 

To obtain approximate values of St as a function of Re 0 and Mq 
for t w /t aw other than 1.0 and 0.5, it is recommended that the results 
of figure 17 he interpolated or extrapolated to give the value of 2St/Cf. 
at the desired condition. Then can he found from equations (34b) 

and (34c) as previously described, and thus the value of St is obtained. 


Relation between Boundary-Layer Thickness and Distance along Plate 


From, the results given in the preceding sections, the skin friction 
or heat transfer for a given boundary-layer or momentum thickness can be 
calculated. In order to calculate the relations between thermal or flow 
boundary-layer thickness and distance along the plate, the well-known 
integral momentum and energy equations may be used. These equations may 
be written as follows for a flat plate (zero pressure gradient): 



(40) 


(41) 


If 6 = 5^ = 0 for x = 0 and equations (40) and (4l) are integrated 
with respect to x, they became, in dimens ionless form, 


Re x 




_P_ 

p w 


u + (uj - u + )dy + 



(42) 


(43) 


where the bracket for the upper limit of integration refers to the value 
of the variable of integration at that point. These equations give the 
relations between 8 + and Re x and 5+ and Re x . The property ratios 

are obtained from equations (10) to (13). 
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Equation (42) can "be written in the more convenient form 



The Reynolds number based on momentum thickness Reg is plotted against 

Re x , as found from equation (44), for an insulated plate in figure 18. 

The value of Reg decreases at a given Re x as the Mach number increases 

if the free -stream properties remain constant. This is caused (eq. (44)) 
by the decrease of friction factor with increasing Mach number (fig. 7). 
Data included in figure 18 agree reasonably well with the analytical 
curves . 

Figure 19 is similar to figure 18, except that t w /t aw = 0.5. For 
given values of Re x and Mach number the values of Reg are generally 
a little higher for t w /t aw of 0.5 than for t w /t aw . of 1. This trend 
can be understood from examination of equation (44), since is higher 

for t w /t aw of 0.5 (fig. 10) than for t^/tg^ of 1.0 (fig. 7). 

Predicted skin-friction coefficients for an insulated plate are 
plotted against Re^j. in figure 20. The trends with Mach number are 
similar to those obtained when Cf is plotted against Reg but are 

less pronounced, because the boundary- layer thickness at a given x de- 
creases with increasing Mach number. Experimental datq, for low-speed 
flow included in the figure are in good agreement with the predicted 
curve for a Mach number of zero. Data for higher Mach numbers are also 
in reasonable agreement with the predicted curves but are seme what more 
scattered than the data in figure 7, where C^. is plotted against Reg. 
This scatter is apparently caused by uncertainty as to the point at which 
the boundary layer actually starts in a supersonic flow. 

In figure 21, the theoretical curves are replotted as Cf/Cf^i 
against Mach number for various Reynolds numbers based on x. The 
effect on G^/C^ i of varying Re x becomes appreciable at high Mach 
numbers . 

Stanton numbers for an insulated plate are plotted against Re x 

for a Mach number of zero in figure 22 . Curves for higher Mach numbers 
and for t w /t aw of 0.5 involve considerably more calculation and were 
not obtained. 

Figures 23 and 24 are analogous to figures 20 and 21, respectively, 
except that they are for t v /t aw of 0.5. The friction factors, as ex- 
pected, are higher for the larger rates of cooling. 
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Reference Temperature for Re x Results 

The customary use of a reference temperature concept req uir es that 
the Reynolds number dependence of Cf and of St be the same for all 

Mach numbers so that Cf/Cf -j_ and St/St^ should not be functions of 
Reynolds number. Examination of the predicted curves of Cf/Cf ^ and 
St/Stf against Mach number as shown in figures 8, 11, 13, 15, 21, and 
24 shows, however, that C^/Cf ^ and St/St^ are strong functions of 

Reynolds number at the higher Mach numbers. Therefore, the present the- 
ory cannot be represented accurately by one reference temperature valid 
for all Reynolds numbers . 

For purposes of comparison, however, the results obtained by using 
Eckert's suggested reference temperature (ref. 3) are shown in figures 
21 and 24. Agreement with the present theory for Re x = 10 8 is quite 
good. Use of Eckert's reference temperature method to solve for Cf/Cf 
is recommended, then, if large values of Re x (near 10 8 ) are considered. 
In order to solve for Cf, the value of Cf ^ can be taken from, the 
curve for Mg = 0 in figure 20. An approximation (within 5 percent) to 
this case is 

' C f = 0.0292 Re x 0-151 

In order to solve for values of Cf at values of t w ./t aw other t han 

O 

1.0 and 0.5 and for values of Re x other than 10 , the same approximate 
procedure as recommended for Reg as the variable can be employed. For 
this case, however, instead of figures 8 and 11, figures 21 and 24 and 
Eckert's reference temperature should be utilized. 

To obtain an approximate relation between Stanton number and Re x 
the following procedure is recommended: Find the value of Reg corre- 

sponding to the specified Re x by interpolation or extrapolation of fig- 
ures 18 and 19. From this value of Reg find the Reynolds analogy fac- 
tor by similar use of figure 17 for the specified values of t w /t aw and 
Mq. This value of the Reynolds analogy factor and the value of Cf ob- 
tained as shown in the previous paragraph are sufficient to solve for 
the Stanton number for the specified conditions. 


Clos ing Re m arks 

Ko attempt has been made in this analysis to include the effects of 
dissociation, shock waves, radiation, slip flow, or induced pressure 
gradients . 
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A rough estimate of the effect of dissociation may be inferred, as 
pointed out by Eckert (ref. 3), from the theory of laminar boundary 
layers. Thus, both Crown (ref. 20) and Moore (ref. 21) conclude that 
the effect of dissociation on friction factor and heat flux for the lam - 
inar boundary layer will be small if the wall temperature is less than 
the air dissociation temperature, which, even at a pressure of 0.0001 
atmosphere, is above 3000° F . Their calculations were made for Mach num- 
bers up to 20. Where dissociation is appreciable, it is recommended 
that the heat-transfer coefficients presented in this report be inter- 
preted as based an an enthalpy difference instead of a temperature dif- 
ference, Thus, 

h = -A! 

Hw - H aw 


and 


St 


(Sw - W U SPS 


where H is the enthalpy. 


Although in practice there would be a shock wave originating near 
the leading edge of the flat plate for high Mach numbers, the effect on 
temperature and pressure distributions appears too complicated to be 
taken into account. Therefore, constancy of free -stream pressure and 
temperature has been assumed. 


The possibility of encountering slip flow at high Mach number must 
also be considered. According to Eckert (ref. 3) the assumption of a 


continuum is valid as long as the Khudsen number 


12 » 
~ “8 


0.499 Re & 


is less 


than 0.01. To obtain a conservative estimate for the range of condi- 
tions considered in this report, values of Mg and Reg of 20 and 10 1 * * 4 , 
respectively, are used. For these values the Khudsen number is 
0.00298, which is well below the criterion for slip flow. 


SUMMARY GF RESULTS 

The following results ■jpire obtained from the analysis of turbulent 
flow and heat transfer over '.a flat plate at high Mach numbers: 

1. The frictional heating that occurs at high Mach numbers produced 

a flattening of the velocity profile, as does heating the plate by other 

means. Cooling the plate caused the velocity gradients near the outer 

edge of the boundary layer to increase. 
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2. The skin-friction coefficients and Stanton numbers at a given 
Reynolds number decreased as Mach number increased. 

3. The curves for the ratio of friction coefficient to the incom- 
pressible coefficient against Mach number agreed closely with the curves 
for the ratio of Stanton number to incompressible Stanton number against 
Mach number. 

4. Cooling the plate to offset the effects of frictional heating 
increased the friction coefficients and Stanton numbers . 

5. Frictional heating at high Mach numbers produced a thinning of 
the boundary layer at a given position on the plate for the same free- 
stream properties. 

6. The predicted friction coefficients and Stanton numbers agreed 
closely with representative experimental data. 

7. The Reynolds number effect on both friction factor and Stanton 
number increases greatly with increasing Mach number. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, January 17, 1958 
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APPENDIX A 

DERIVATION OB' TURBULENT MOMENTUM AND ENERGY EQUATIONS 
Momentum Equation 

The momentum equation for compressible boundary- layer flow past a 
flat plate can be written as 


du Su b ( du\ 

pu 3E + pT 3? ‘ ^ ( p Sf) 


and the continuity equation as 


+ d(pv). = 0 

ox oy 


(Al) 


(A2) 


Time derivatives and pressure gradients are neglected in equations (Al) 
and (A2 ) , as they drop out when time averages are taken. 

The instantaneous quantities in equation (Al) are now replaced by 
their time averages and fluctuating components, which are written as 


u = u + u' 
v = v + v' 



(A3) 


and time averages are taken term by term. The following order-of- 
magnitude criteria are used for both the momentum and energy equations: 


■£ic 


- 0 ( 1 ) 


|f ■ °C £fl > 

Pi u, t - 0(1) 
v * 0(5) 
k, p - 0(S 2 ) 
u'v ' , p 'v', etc. * 0(6) 

o-f =« o(s 3 / 2 ) 


p’u'v’, p'u' 1 *, etc. 


Double correlations containing k* and p 1 * 0(S ) 
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The first five of these criteria are the usual boundary- layer assumptions. 
The sixth results from assuming that the laminar and turbulent shear 
stresses are of the same order of magnitude. The seventh is consistent 
with the sixth, since a triple correlation should be roughly of the mag- 
nitude of a double correlation raised to the 3/2 power. The eighth 
appears justified since it might be expected that k' and _p ' should 
be at least one-half order of magnitude less than k and p. 

With the preceding criteria, the time -averaged momentum equation 
becomes, on neglecting terms of magnitude 8 and less, 

p u + (57 + 777) |p = |s - 9 777 ) (A4) 

and the continuity equation, 

^ ^ - 0 C*> 

Considering the relation 

pv + p 'v' = pv (A6) 

equations (A4) and (A5) can be rewritten 

55 s + ^ ( 5 H ‘ p ' r7T ) (A7) 

and 

+ mp. = 0 (A8) 

C caparison of equations (Al) and (A7) leads to the definition of t as 

t = p ^ - p u’v’ (A9) 


Energy Equation 


The energy equation 
flat plate is 

&fc 

e uc p37 


for compressible boundary- layer flow past a 

dt 3 f ' St\ /’du\ 2 


(A10) 
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Time derivatives and pressure gradients are again neglected, as time- 
averaging cancels them. If the momentum equation (Al) is multiplied 
through by u and then added to equation (AlO), the result is 


pu 


a 

x 





where Cp 


is considered constant. 


(All) 


Again substituting for the instantaneous quantities the sum of the 
time-averaged and fluctuating components, and then neglecting terms of 
small order of magnitude on taking time averages, equation (All) becomes 

p “ s (v + t) + ^ + P 777 ) (“s* + t) - 


(A12) 


(A13) 


(A14) 

It should be noted that this treatment gives no density fluctuation 
terms in the expressions for t and q. This same result was found by 
Van Driest (ref. 22) and by Rubes in (ref. 19). Combining p'v' with 
pv and writing the_sum as pv present no difficulty, because, in a 
complete solution, pv could be el im inated from the momentum and energy 
equations by the equation of continuity. An assumption for p'v' would 
be necessary only if it Were desired to calculate v. 


9 |r 3 1 , — Su — — TT-T 

- PC P V * 




pu 


u'v’^ 


Again employing equation (A6), equation (A12) becomes 


-- a 
pu S 


(pp* + t) 


+ p7 If ( V + v) 


3 fr-bt . - - Su — 


#(' k $ + - pc p v,t ' - pu 

A ccmpariscn of equations (All) and (A13) shows that 


u’v'^ 


-( 


r- St — Su — — rrr — 

+ “ pc p v * - pu 


u'v^ 
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APPENDIX B 


PROOF THAT — = JL P® Pr = a = 1 


w 


*w 


In terms of eddy dif f us ivit ies , the momentum and energy equations 
may be written as follows: 


— - du , — c)u 5 (— - > 5u 

pu ^ + pv ^ = -S5F pe ^ 


-- a 

pu 


s (v + t) + ^ b + 4 ) = 


3y 


& + P®p«h^ + uGT + pe) ^ 

(B2) 


(Bl) 


— _\ c)u 


The energy equation (B2) can be rearranged to read 


u 

7 


— 5 
pu ^E 


(v + t) + pv If (v + x) 


= _p? + P £ h ^ Cc p t) + 01 + pe) ^ ^ 


(BS) 


For Pr = a = 1, equation (B3) becomes 


p{1 azr " 


(v + f ) + ^ ^(v + t) - ^ <5 + + vj < E4 > 


If equations (Bl) and (B4) are each solved for pu and the results 
equated, there is obtained 


B r- , — v Su -- du 3 
|_g 1- pe) - pv ^ ^ 




C5 + + t) -JZ. + t! 


Su 

!& 


b 


(v + v) 


The assumption is now made that 


_ -2 _ 

Cpt + i- = Au + B 


(B5) 


(B6) 


_ -2 

The expression given by equation (B6) for c^t + is substituted into 

the right side of equation (B5) . Since the right side becomes identical 
with the left side upon this substitution, equation (B6) is a valid 
relation . 
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The constants A and B In equation (B6) are evaluated as follows: 
at y = 0: 

u - 0, t = B = ^ 

at y S3 0: 

77 du _ _ r 5t A %• 

^ " v ' k A = - -r 

With the constants thus evaluated, equation (B6) becomes (dropping the 
bars for convenience) 

u 2 %r 

C t + _ 5= - — u + c t , 

P 2 x w p w (B7 ; 

If equation (B7) is made dimensionless, it becomes simply 


u + = & 


(B8) 


Previously obtained relations for t/t w and q/q w are 





JL - fJL 

. P- 1 

\ du + 


(5) 





Pw Pw/Pw^ 

/ =3y+ 


_S_ « 


1 

+ JL a € , ^ 

dt + a 

u + ( JL 4 . J2_ e , \ 

du + 

( C\ 

%r 

K 

Pr 

w 

Pw Pv/PwJ 

dy+ ’ P 

\Pv Pv Pw/PvJ 

dy+ 


If use is 

made 

of equations (B8) ; 

, (lO), and (26) and the fact 

that 



S. _ JL 


(B9) 


QTQ* 
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APPENDIX C 


VELOCITY MD TEMPERATURE HtQF ILES FOR CONSTANT RATIO GF TURBULENT 
TO MOLECULAR SHEAR SffiESS AT yj AND MOLECULAR SHEAR S1RESS 
AND HEAT TRANSFER IN REGION AWAY FROM WALL (Pr = l) 


The equation for velocity profile used near the wall is equation 
(16), where 1 - at + ' = 1 - au+2 for Pr = 1 and 6 = 0 (see eqs. (B8), 
(26 ) ) and definitions of t + and t+'). 


Previously, the expression for t/ t w was shown to he 

jr ^ du + _p_ e du* 

T w dy + p w ^wT^w dy + 

and far from the wall (eq. (8)), 



( 5 ) 


Combining equations (5) and (8) and assuming constant shear stress across 
the Boundary layer give 



(Cl) 


The variations of density and viscosity with temperature are 



J 


(C2) 
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For Pr = a = 1 and P = 0, the temperature ratio is expressed as 


qr- = 1 - ccu +2 


Substitution of equations (C2) and (C3) into (Cl) yields 

,+\3 

„ (a 

1 = I (l - au +2 ) 


)0.68 + x 2 U + i 
1 - tm +2 / d 2 u+ y 

W +2 J 


du 


dy* 


Solving for d 2 u + /dy +2 gives 


2,i+ 


d^u 


-W 


<3y 


.+2 


\/ (1 - 0.1+ 2 ) |l - (! . cca+2) 0 - 68 1 >£ 

V [_ dy+ 


If a change in variables is made as 


v = 


du 

dy + 


equation (C5) can be integrated to give 

vU+ 


/: 


xdu* 


v = v 1 e 


/^(1-cgu +2 ) [I- (l-au +2 )° * 68 v] 


(C3) 


(C4) 


(C5) 


(C6) 


j_ -4- 

The solutions for u as a function of y can be obtained by a process 
of iteration. Assumed values of v for a given incr eme nt in u + are 
substituted into the right side of equation (C6) until it eq ual s the 
left side. The relation between u + and is then calculated from 



(C7) 


Fran equation (B8), u + = T" 1- for q/q w = t/t w , so that the relation be- 
tween T* and y + is also known. 


OT^ • 
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APPENDIX D 


VELOCITY PROFILES FOR LINEAR VARIATION OF SHEAR STRESS 
AND HEAT TRANSFER ACROSS BOUNDARY LAYER (Pr = l) 

Because the variation of shear in the thin region near the wall is 
negligible, the same equations are used in the present case as were used 
in appendix C. 

From equation (5), the equation for t/t v far from, the wall, neglect- 
ing the viscous stress, is 


_r g_ e du + 

T w p w *V p w dy + 


The expression for e/ (p^/p^.) far from the wall is, from equation (8), 



For a linear variation in shear stress. 



(Dl) 


Combining the foregoing equations yields 



For p = 0, Pr = a = 1 (t/t^ = q/%), 



(D2) 


(D3) 
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Substituting equation (D3) into (D2 ) , rearranging, and taking the square 
root of both sides give 


d 2 u + _ - x 

( 3 y +2 ^1~". ou+2 



(D4) 


Letting v = du + /dy + 


in equation (D4) and integrating give 



(D5) 


Equation (D5) can be solved by iteration for [3 = 0 and a given a and 

S + , to give u+ as a function of y + . From equation (B8), u + = T* for 
q/q w = t/t so that the relation between T+ and y + is also known. 
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Figure l. - Generalized computed velocity distribution for constant-property turbulent flow (ref. 4) with Constanta in 
equations determined from pipe data conpared with data for low-speed boundary-layer flow on a flat plate. 
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Figure 2. - Predicted effect of variation of shear stress and heat transfer across bo und ary layer on velocity and tempera- 
ture distributions. Prandtl number, 1; heat-flux parameter P, 0. 
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Figure 3. - Effect of various assumptions for variation of yj on generalized velocity or temperature distribution. 
Prandtl number, 1; heat-flui pars peter p, 0. 
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Figure *. - Predicted generalised velocity distribution for air vith heat transfer and frictional heating. Prandtl if». 

number, 0.73. H 
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Figure 5. - Predicted generalized total -te up erature distribution for air vith heat transfer and frictional beating. 
1 Praadtl nuaber, 0.73. 
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Figure 7 , - Tariation of predicted sfcin- friction coefficient vith Btttneptm-thiclmeia Reynolds umber *nfl Mach 
umber for insulated plate and copper! son vith experiment* Prandtl number, 0.73; heat-flux parameter 3, 0- 
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Figure 8. - Variation of Cf/Cf ^ with Mach number for various values of momentum-thickness Reynolds number for insu- 
lated plate. Prandtl number, 0.73; heat-flux parameter P, 0. 
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Figure 9. - Variation of Cf/Cf ^ with Mach number for momentum-thickness Reynolds 
number of 6000 for insulated plate and comparison with experiment . Prandtl humber 
0.73; heat-flux parameter g, 0. 
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Figure 11. - Variation, of C f /C fjl 
ty/taw of 0.5. Prandtl Timber, 


vith tfcoh number tar tbjIoue values of ucwentuo-thlctoeBB Reynolds mmtoer and 
0.73. 
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Figure 13. - Variation of St/sti vith Mach number for various values of ocnuentuoi-tblclmese Reynolds number for insu- 
lated plate. Erandtl number, 0.73; beat-flux parameter 0, 0. 
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Figure IS. 
Of 0.5. 


- Variation of St/Sfc^ with Mach rraafoer for various values of rocmentf -t h ic K neas Reynold# lumber end t„/t av 
Praudtl nimher, 0.75. 
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Figuxe 16. - Predicted variation of temperature -recovery factor with momentum-thickness 
Reynolds number and Mach number for e^/e =* a = 1 -and 1.07. Prandtl number, 0-73- 
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Longitudinal-distance Reynolds number, Re x 

Figure 18. - Predicted variation of momentum-thickness Reynolds number with longitudinal- 
distance Reynolds number mi Mach number for insulated plate and comparison with experi- 
ment. Prandtl number, 0.73; heat-flux parameter (3, 0. 
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Figure 19. - Predicted variation of momentum-thickness Reynolds number with longitudinal- 
distance Reynolds number and Mach number for t^/t^ of 0.5. Prandtl number, 0.73. 
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Figure 20. - Predicted variation of akin- friction coefficient with lcngitudinal-distance Reynolds number and Mach number 
for insulated plate and coejparlaon with experiment . Prandtl number, 0.73; heat-flux parameter p, 0. 
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Figure 22. - Predicted variation of Stanton number with longitudinal-distance Reynolds num- 
ber for insulated plate and comparison with experiment for low-speed flew. Prandtl num- 
ber, 0.73J Mach number, 0. 
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Figure 25. - Predicted variation of skin- friction coefficient with longitudinal-distance Reynolds number and Mach, number 
for tw/tav of 0.5. Fraud tl number, 0.73. 
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Figure 24. - Variation of Cf/Cf^i with Jfech maabar for various values of longttuiinal-Aistanas Keyaolds mtmbar and 
tw/' t e.w o:f 0- 5 - Prandtl rrunber, 0-73. 
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